Barkhausen avalanches in anisotropic ferromagnets with 180° domain walls 



Bosiljka Tadic 1 '* and Ulrich Nowak 2 -** 
1 Jozef Stefan Institute, P.O. Box 3000, 1001 Ljubljana, Slovenia 2 Theoretische Tieftemperaturphysik, 
Gerhard-Mercator-Universitat-Duisburg, D-47048 Duisburg, Germany 

We show that Barkhausen noise in two-dimensional disordered ferromagnets with extended do- 
main walls is characterized by the avalanche size exponent r s = 1.54 at low disorder. With increasing 
disorder the characteristic domain size is reduced relative to the system size due to nucleation of 
new domains and a dynamic phase transition occurs to the scaling behavior with r s = 1.30. The 
exponents decrease at finite driving rate. The results agree with recently observed behavior in 
amorphous Metglas and Fe-Co-B ribbons when the applied anisotropic stress is varied. 



Barkhausen noise measured in disordered ferromagnets 
at low temperatures under the condition of slow driv- 
ing by an external field through the hysteresis loop ex- 
hibits scaling behavior without tuning of any parameter 
The scaling properties of Barkhausen avalanches 
obtained in various alloys can be grouped in three dis- 
tinct universality classes differentiated by the value of 
the avalanche size exponent as § r s « 1.3, r s « 1.5, 
and t s « 1.7. However, the origin of the scaling of 
Barkhausen noise (BN) and the occurrence of different 
universality classes is still not fully understood. 

The above mentioned measurements jlj-ffl] are done 
mostly on thin ribbon samples (thickness d < 200/^m). 
The variety of measured scaling exponents can be at- 
tributed to the differences in the applied driving condi- 
tions and to the diversity of domain structures occurring 
in different samples. It was found that the exponents de- 
crease continuously with increasing driving rate [^,[5],^|. 
The structure of domains in commercial alloys has been 
studied by variety of techniques (see for instance |9]). In 
thin ribbons of the amorphous Metglas Fe-B-Si [l|,0j9| 
and Fe-Co-B alloys |||| which are annealed in a paral- 
lel field H or under an applied anisotropic stress ||[7]] 
a structure with few domains occurs with 180° domain 
walls parallel to the anisotropy axis. The demagnetiz- 
ing fields, which depend on the form of the sample, and 
the range of interactions play an important role for the 
equilibrium domain structure as well as for the domain- 
wall dynamics |5H|. It has been demonstrated that the 
domain structure in the amorphous Metglass J?]] and Fe- 
Co-B alloys || can be controlled by varying tensile stress 
and that short-range interactions dominate over dipolar 
forces Q. In addition, in these systems the demagne- 
tizing fields are minimized with longitudinal anisotropy 
||. Hence, in our study we can neglect dipolar forces, 
but we take into account the existence of domain walls. 
Numerical studies of BN using short-range Ising models 
with various types of disorder ||lo|-|i"2| (see also Jl3|] for 
more realistic type of interactions) usually start from a 
uniform ground state and nucleate a random pattern of 
clusters of reversed spins by increasing the external field, 
thus neglecting the preexisting domain structure. 

In this paper we simulate Barkhausen avalanches us- 



ing a model with preexisting extended domain walls which 
are confined in two dimensions, motivated by the stress- 
induced anisotropy in realistic systems [HJtJ. Increasing 
the external magnetic field the domain wall may either 
move through a random medium or new domains may 
nucleate when this is energetically favorable fiij . It has 
been suggested recently that interface depinning in a ran- 
dom medium is responsible for the scaling behavior of BN 
however, models of an elastic interface yield scal- 
ing exponents which are lower than the measured ones 
PJl- In this work we use a ferromagnetic model with 
short-range interactions and a random-field pinning and 
we show that the results compare well with two univer- 
sality classes of measured scaling exponents. In addi- 
tion, a dynamic phase transition between these two scal- 
ing behaviors appears when the strength of disorder (or 
size of the domains) is tuned. It has been recognized 
recently that disorder effects are enhanced leading to 
smaller domain sizes by decreasing either tensile stress 
HQ or grain sizes ||l3). Motivated by these suggestions, 
here we study the influence of disorder on the scaling 
properties of Barkhausen noise. We expect that the re- 
sults are relevant to realistic samples in which the domain 
size exceeds the sample thickness. 

We consider an Ising model on a square lattice of size 
L x L assuming that local random fields hi are generated 
by coarse-graining from an original disorder ||14| in the 
presence of the external magnetic field H: 



U = - J i,jSiSj - ^2(hi + H)Si 



(1) 



<»,i> 



Here we set Jjj = 1 to be a constant interaction be- 
tween nearest-neighbor spins Si = ±1. Hence, all fields 
and energies are measured in units of Jij. A Gaussian 
distribution of hi is assumed with zero mean and width 
/. We create an initial domain wall in (11) direction 
by a rotation of the lattice by 7r/4 and setting all spins 
except of those in the first row opposite to the exter- 
nal field [jl5| . Periodic boundaries in the direction of the 
interface and fixed boundaries in the perpendicular direc- 
tion are applied. The dynamics consists of a spin align- 
ment parallel to the external field when the local field 



1 



h!f c = V ■ JijSj + H + hi exceeds zero. For the simula- 
tion of a hysteresis loop the field updates are adjusted to 
the minimum local field (infinitely slow driving), thus the 
driving field is uniform in space but fluctuates in time. 
We also briefly discuss the effects of finite driving rates. 

It should be stressed that at each time step we update 
all spins, which is suitable for a globally driven magnetic 
system. In this way new domains are nucleated when it 
is energetically allowed, in contrast to models of driven 
interfaces, where an update is restricted only to the sites 
located next to the interface, resulting in "percolation" 
growth even at high disorder [ jl6| . Another important 
feature of our model is the anisotropy |l7j due to the ini- 
tial conditions: the threshold driving forces in the parallel 
and perpendicular directions appear to be different. This 
model has the following advantages: (1) In the limit of 
vanishing disorder the (11) interface can be moved by an 
infinitesimally small field. This bypasses the problem of 
threshold energy 2J, which occurs in an (10) interface 
p6 | implying that a large lattice size have to be used in 
order to find a spin with a random field large enough to 
surmount the energy barrier. Therefore, here we can ap- 
ply smaller lattice sizes and vary disorder configurations 
(we use up to 2 x 10 3 configurations and up to L = 400). 
(2) The (11) interface depinns at infinitesimally small 
field at low disorder / — » 0. Hence, an upper limit of dis- 
order f* exists at which depinning is no longer possible, 
and nucleation of new domains in the interior becomes fa- 
vorable at large enough fields. This feature comprises an 
important difference compared to the model of Ref. jlO) , 
where nucleation of a single spanning cluster at / — > 
requires energy threshold 4 J, and thus becomes obscured 
by finite lattice size and large fields. 

Fig. 1 shows snapshots of simulated systems for high 
and low disorder. For low disorder only domain wall mo- 
tion occurs. For values of / that exceed a certain critical 
value /* (determined below), domains nucleate inside the 
system, thus leading to the same structure of clusters as 
in systems without an initial domain wall Jl2|,[l9| . These 
two regions are shown in the phase diagram in Fig. 2 
together with simulation results for the coercive fields 
Ho(f) of the hysteresis and the critical fields H c (f) of 
the depinning transition (explained below). 

We apply the quasi-static driving described above (an 
example of time series of field increments is given in the 
inset to Fig. 2) and monitor the motion of domain walls. 
The number of flipped spins s between two consecutive 
locally stable configurations of the wall determines the 
size of Barkhausen avalanche. The avalanche size dis- 
tribution D(s, f) is shown in Fig. 3 for various values 
of disorder /. Taking the avalanche statistics along the 
ascending part of the hysteresis loop (until eventually de- 
pinning occurs) and averaging over many disorder con- 
figurations we find the slopes according to D(s) ~ s~ Ts 
as t s = 1.54 for / < 0.6, and t s = 1.30 above /* w 0.6. 
The cut-off also decreases with /. A similar behavior 
was found experimentally in Metglas 2605TCA in Refs. 
M . With the stress a varying in the range from - 525 



MPa the slopes of the size distributions are reported to 
vary from 1.29 to 1.60, and the cut-offs increase 0. In 
Fe64Co2iBi 5 the measured slope was 1.28 Q for a up to 
140 MPa. The applied anisotropic stress stretches the 
domain walls, thus for the degree of disorder / we have 
/ ~ l/a x . We also measure the distributions of the linear 
extensions of avalanches in the directions parallel (w) and 
perpendicular (h) to the wall leading to the anisotropy 
exponent ( = (t w - l)/(r, l - 1) = 0.92. (See Table 1). 




FIG. 1. Emerging cluster structure of flipped (bright) and 
unflipped (dark) spins, and most recent cluster (gray), illus- 
trating self-afnne growth at low disorder (bottom), and nucle- 
ation and blocking by previous clusters at high disorder (top) 
L = 200. 

In the low disorder region the domain wall depinns 
when the driving field exceeds a critical value H c (f). To 
investigate the depinning transition we start with a flat 
wall and apply a constant field H measuring the long- 
time limit of the domain wall velocity u, which is defined 
as the number of flipped spins per internal time step rel- 
ative to L. The velocity averaged over 2000 disorder 
realizations, [v], is the order parameter of the depinning 
transition. Here [v] and its fluctuations are obtained for 
L =100, 200, and 400 and the critical exponents and crit- 
ical fields H c (f) (also shown in Fig. 2) are determined by 
a finite size scaling plot of the general form 



Y(X,L) = L- ZY y(L 1/u X) 



(2) 



with X = (H — H c (f))/H c (f) (see Table 1, left side, and 
Ref. JnJ). Here Zy = [3/v for the analysis of the order 
parameter [v], and Zy = j/v for the analysis of the fluc- 
tuations of [v] . Our results are compatible with the corre- 
lation length exponent v = 1.23 ± 0.04, (3 = 0.43 ± 0.03, 
and 7 = 1.52 ± 0.06. The exponents are universal in 
the region 0.35 < / < 0.6. Using the scaling relation 
[3/v = z — Q valid for the depinning transition Jl7t we 
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find the dynamic exponent z =1.27. Then the relation 
D(t s — 1) = z(r t — 1), where the fractal dimension of 
anisotropic avalanches is D = 1 + £, leads to the du- 
ration exponent r t =1.83. These results are in good 
agreement with universal criticality in a class of driven 
dynamical systems recently discussed in Ref. The 
fraction of active sites at time t scales as N(t) ~ t K , 
with k = D/z — 1 = 0.51, compared with 0.58 in 
The value for the exponent v compares well with one in 
|P^ , ^8| . For the (local) roughness exponent ( values in 
the literature ranging from 0.5 to 1.23 can be attributed 
to the influence of the anisotropy [jl7| , elasticity of the in- 
terface [[22), distribution of disorder |L5 16 1, and driving 
conditions p0j , ^3| . Further analysis, e.g., by the dynamic 
renormalization group, is necessary in order to determine 
if the value £ = 0.92 found in this work represents a new 
universal behavior or a crossover due to finite size effects. 
This value suggests closeness to the class of models with 
vanishing interface velocity |is}| and self-organized depin- 
ning p3| . At low disorder lateral motion of the interface 
makes it possible to overcome strong pinning centers and 
to maintain the self-affine growth (see Fig. 1). However, 
at the transition we find cj> = 2j/(1 — £) — (3 = —0.23, 
thus A e // — y resulting in a different critical behavior 
compared to the case A e // — > oo studied in Refs. JTt| , |18| ] . 

In order to determine the largest disorder where de- 
pinning is still possible, /*, we can formally extend the 
above study of the domain wall velocity now averaged 
over hysteresis loop and disorder [< v >], and the cor- 
responding susceptibility. Applying then Eq. (|^) with 
X = (/ — /*)//* (scaling collapse is shown in the inset 
to Fig. 3) we find /* = 0.61 ± 0.02 and the exponents /3, 
7 and v shown in Table 1, right side. 

In the region of high disorder / > /* nucleation of 
new domains of finite size which are blocking each other's 
spatial extent becomes the dominant feature which deter- 
mines the scaling properties of BN (see Table 1 (right)). 
The section of power-law behavior of the distribution 
of avalanches increases with decreasing /, in qualitative 
agreement with increasing stress in experiments |^,|tJ . In 
the inset to Fig. 4 the distribution of avalanche durations 
is shown for various values of / above /*. The main Fig. 
4 shows the scaling plot P(tJ) = (Sf) ZVTt V(t{Sf) zv ), 
where Sf = (/ — /*)//*, obtained by using /* = 0.62 and 
the products of the exponents zvr t = 4.43 and zv = 2.98. 
^From a similar scaling collapse of the size distribution 
we find Dvt s = 5.61 and Dv = 4.30. Thus, these val- 
ues together with the ones in Table 1 (right) lead to 
v = 2.3 ± 0.1, 13 = 0.12 ± 0.03, and 7 = 5.2 ± 0.1. 
Note that the value of /* = 0.62 within statistical er- 
ror bars ±0.03, which we estimated from two different 
types of scaling fits (see Figs. 3 and 4) is by no means 
definitive. (The error bars are expected to increase when 
more scaled quantities or wider range of system sizes are 
explored). However, since no exact results are available, 
this value can be regarded as a rough estimate of the 
critical disorder. 

By applying a finite driving rate r = AH/ H maxi where 



is the saturation field, the cut-offs of the distribu- 



tions in the high disorder region increase (see also (L2j). 
Whereas slopes of the distribution decrease with r due 
to mainly the coalescence of avalanches. For the size of 
avalanches we find, for instance for / = 0.88, r s = 1.18 
for r = 0.01, and r s = 1.10 for r = 0.02. 

In conclusion, we have shown that the anisotropic 2- 
dimensional motion of domain walls pinned by quenched 
impurities and short-range interactions are relevant for 
the scaling behavior of Barkhausen noise in the pres- 
ence of extended domain walls. We find two universality 
classes with t s = 1.54 and t s = 1.30, in a good agreement 
with experiments in amorphous Fe-B-Si and Fe-Co-B rib- 
bons under anisotropic stress. By slow driving the scaling 
behavior is robust in a wide range of the driving field and 
disorder (or stress) values. These two universality classes 
correspond to the motion of extended domain walls (i.e., 
at low disorder or high stress), and many finite domains 
(high disorder or low stress), respectively. The scaling 
exponents decrease with finite driving rate. We find that 
the domain structure changes via a dynamic phase tran- 
sition at a critical disorder, which can be directly moni- 
tored in experiments by tuning uniaxial stress. 
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Exponent 


Value (LD) 


Exponent 


Value (HD) 


T s 


1.54 




1.30 




2.03 


n 


1.47 


Th 


2.12 


D 


1.88 


c 


0.92 


z 


1.23 


P/v 


0.35 


13/u 


0.06 


j/v 


1.22 


j/v 


2.26 


1/v 


0.82 


l/u 


0.42 



TABLE I. Scaling exponents obtained by direct simulation 
of various quantities for the low (LD) and high (HD) disorder 
region. Error bars estimated from independent fits are within 
±0.03. 
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FIG. 2. Phase diagram in the plain disorder-driving-field 
(both in units of J). The vertical dashed line at /* = 0.62 sep- 
arates high and low-disorder region. Long-dashed line with 
crosses: coercive fields Ho{f); filled symbols: critical fields 
H c (f) obtained from scaling fits. Inset: Time (in MCS) se- 
ries of driving field increments AH (in J/MCS) for / = 1.2. 
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FIG. 3. Probability distribution of avalanche size D(s, f) 
vs size s for various values of disorder /=0.4, 0.6, 0.9, 1.0, and 
1.1 (bottom to top). Inset: Scaling collapse according to Eq. 
(0) of the velocity of domain wall [< v >] vs. / for L=128, 
196, and 256. Each point is averaged over 200 samples. 
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FIG. 4. Inset: Probability distribution of durations of 
Barkhausen avalanches P(t, f) vs. duration t in MCS, for 
L = 400 and for /= 1.2, 1.0, 0.96, 0.92, 0.88, and 0.84 (left to 
right). Main figure: Scaling collapse of the data (see text). 
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